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Motivated by the problem of Many-Body Localization and the recent numerical results for the level 
and eigenfunction statistics on the random regular graphs, a generalization of the Rosenzweig-Porter 
random matrix model is suggested that possesses two transitions. One of them is the Anderson 
localization transition from the localized to the extended states. The other one is the ergodic 
transition from the extended non-ergodic (multifractal) states to the extended ergodic states. We 
confirm the existence of both transitions by computing the two-level spectral correlation function, 
the spectrum of multifractality f{a) and the wave function overlap which consistently demonstrate 
these two transitions. 


I. INTRODUCTION 

Motivated by the problem of Many-Body (MB) 
Localizatiorf^ and the applicability of the Boltzmann’s 
statistics in interacting disordered medicP, there was re¬ 
cently a revival of interest to the Anderson localization 
(AL) problem on hierarchical lattices such as the Bethe 
lattice (BL) or the random regular graph (RRG). Due to 
hierarchical structure of the Fock space connected by the 
two-body interaction, statistics of random wave functions 
in such models is an important playground for MB local¬ 
ization. In particular, the non-ergodic extended phase on 
disordered hierarchical lattices could model a breakdown 
of conventional Boltzmann statistics in interacting MB 
systems and an emergence of a phase of a “bad metal”!^ 
or unconventional fluid phased in systems of interacting 
particles. 

However, even for the one-particle AL existence of such 
a phase in a finite interval of disorder strengths is a highly 
non-trivial issue. 

According to earlier studie^^ there is only one tran¬ 
sition in such models at a disorder strength W = Wat 
which is the AL transition that separates the localized 
and ergodic extended states. However, recent numer¬ 
ical studied of level statistics on RRG seem to indi¬ 
cate on the second transition a.t W = Wet < Wat 
which is identified as the transition between the ergo dic 
and non-ergodic extended states. Subsequent studie^^ 
raise doubts about the existence of the second transi¬ 
tion on RRG. Numerical results of Ref.indicate on the 
non-ergodic states on RRG in a wide range of disorder 
strengths down to very low disorder W = 5 -€1 Wat ~ 
17.5, while in Ref.^^ it is demonstrated how an apparent 
non-ergodic behavior for the intermediate matrix sizes 
N in Levy Random Matrix (RM) ensemble evolves into 
the ergodic one at larger N's. Gomplexity of RRG and 
the controversy associated with existence of the ergodic 
transition a.tW = Wet necessitate a search for a simpler 
model in which such a transition may occur. 

Inspired by the success of Wigner-Dyson RM theorjIiSl 


which predictions are relevant in such seemingly differ¬ 
ent fields of physics as nuclear physics and nano- and 
mesoscipic physics, our goal is to search for a RM model 
that would be able to give a simple and universal de¬ 
scription of all the three phases: good metal, MB insula¬ 
tor and “bad metal”, which are relevant in the problem 
of MB localization. An important heuristic argument to 
construct such a model is that RRG with disordered on¬ 
site energies £i is essentially a two-step disorder ensemble. 
The disorder of the first level is the structural disorder 
due to the random structure of RRG where each of N 
sites of the graph is connected with the fixed number 
AT -b 1 of other sites in a random manner. An ensemble of 
tight-binding models on such graphs with deterministic 
on-site energies £i and hopping integrals is believed to be 
equivalent to the Gaussian RM ensembl^^il. The disor¬ 
der of the second level is produced by randomization of £i 
fluctuating independently around zero with the distribu¬ 
tion function p{e). For numerical calculations this distri¬ 
bution is often taken in the form p(e) = W . q (W. _ 
with 9{x) being the Heaviside step function. 

One can expect that the follo wing RM ensemble (the 
Rosenzweig-Porter (RP) ensemblei^ is a close relative of 
RRG with on-site energy Si disorder. It is an ensemble 
of N X N random Hermitian matrices which entries Hnm 
with n > m are independent random Gaussian num¬ 
bers, real for orthogonal RP model (/3 = 1) and complex 
for the unitary RP model (/3 = 2), fluctuating around 
zero with the variance (|i7nmP) = (,3/2) cr. The diago¬ 
nal elements have the same properties with the variance 
i^nn) = 1- The case cr = 1 corresponds to the Gaus¬ 
sian Orthogonal (GOE) or Gaussian Unitary (GUE) en¬ 
sembles and represents the structural disorder in RRG. 
The additional Ci-disorder in RRG corresponds to tr < 1. 
One should also take into account that in order to sig¬ 
nificantly deviate from the GOE or GUE behavior, the 
ratio (|i7„mP)/(77^„) must be proportional to some neg¬ 
ative power of the matrix size N, as the number of the 
off-diagonal terms is ~ A^ times larger. Thus we consider 
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the model: 

(i?!) = 1, = {P/2) a = \^/N\ (1) 

where A ~ 0{N^) and 7 is the main control parameter 
of the problem. One can estimate the strength of dis¬ 
order required for the Anderson localization transition 
as corresponding to the typical fluctuation of diagonal 
matrix element equal to the typical off-diagonal matrix 
element times the coordination number K. For the coor¬ 
dination number K ^ N (each site is connected with any 
other one) this results in y/aN ^ 1, or cr ^ l/N"^. How¬ 
ever, this estimation does not take into account a ran¬ 
dom, sign-alternating character of the off-diagonal matrix 
elements. It is likely that for sign-alternating hopping 
there is another relevant coordination number ^ y/N 
with the critical scaling a ^ 1/N. As we show below 
it corresponds to the ergodic transition. For technical 
reasons the most significant prog ress in the analytical 
studies of the model was achieved ^^ l ^"^ * ^^ ! for the “uni¬ 
tary” RP (URP) ensemble. The conclusion was that at 
7 = 2 the spectral form-factor (two-level correlation func¬ 
tion) isneither of the Wigner-Dyson nor of the Poisson 
forn j^^ l ^^ -11^ which is typical for the AL transition point. 
In contrast, at 7 = I t he leve l statistics was found to 
be GUE!S1. The paper^i^E^E^ have a status of classic 
keynote papers in the field. 

The Dyson ideas of the Brownian motion of energy 
levels first applied to the RP ense mble in Ref.E^were de¬ 
veloped in the series of work^^^EH, Jt was shown that the 
possible transitions in the level statistics are associated 
with the fixed points of parameter A = a{N)/[6{N)]'^ oc 
N~'^/[6{N)]'^, where S{N) is the mean level spacing. 
Then assuming S{N) (x 1/N established in Ref. ESI for 
7 > 1 one obtains the transition point at 7 = 2. If, 
however, the Wigner-Dyson semicircle level density is as¬ 
sumed with S{N) oc 1/VN, then the transition would 
occur at 7 = lE3. Unfortunately, S(N) oc I/Vlv only at 
7 = 0. For 0 < j ^ 1 the following result is valid (see e.g. 
Eq. (170) in Ref.ESl) for the mean density of states p(A) = 
V2S-Ey{TrS), where 5 = E„{\Hnm?) oc Thus 

we obtain S{N) = l/{Np{0)) ex A^“(i+t')/ 2 for 0 < 7 < 1 
and S{N) ex 1/A^ for 7 > 1 , resulting in A oc for 
0 < 7 < 1 and A oc N'^~o for 7 > 1. We conclude that 
the only fixed point of A is possible at 7 = 2, and no 
transition at 7 = 1 can be obtained from the results of 
Refs.ElEll. 

In this paper by a more sophisticated analysis of the 
two-level correlations and the eigenfunctions statistics we 
show that the above extension of the Rosenzweig-Porter 
model indeed contains not one but two transitions. One 
of them at 7 = 2 corresponds to the transition from the 
extended to the localized states. However, the extended 
states emerging at 7 < 2 are not ergodic: their support 
set contains infinitely many sites in the N ^ 00 

limit, which, however, is a zero fraction of all sites, since 
Di < 1. Such non-ergodic extended states on RRG are 
recently discussed in Ref.® With further decrease of 7 


the second transition at 7 = 1 happens which is a transi¬ 
tion from the non-ergodic extended states to the ergodic 
extended states with Di = 1 similar to the eigenstates of 
the GOE. 

We prove this statement in three steps. As the first 
step we use the perturbative arguments to compute the 
statistics of wave function amplitude \ip{ro)\'^ in a cer¬ 
tain observation point Tq. We obtain a drastic change 
of the character of this distribution at 7 = 1 and 7 = 2 
which is summarized in Fig. This result is fully con¬ 
firmed by a numerical diagonalization of the Hamiltonian 
(see Figs. It is also confirmed by the numerical 

analysis of the moments of random wave functions which 
determine their Shannon entropy and the support set di¬ 
mension Di (see Fig. [^. Then we compute numerically 
the overlap of amplitudes for two different wave functions 
with the energy difference w and find the scaling with N 
of the Thouless energy Exh ~ which exponent z 

changes abruptly at 7 = 1 and 7 = 2 . Finally, we per¬ 
form a rigorous calculation of the spectral form-factor 
which also shows the transition at 7 = 1 and 7 = 2 (see 
Fig.0. In the last section we compare the corresponding 
results for our model and for the RRG and demonstrate 
their similarity. It allows us to unify both models in a 
special universality class of random hierarchical models 
which differs from the one realized in localization transi¬ 
tion points of two- and three-dimensional Anderson mod¬ 
els. Further details concerning this model can be found 
in Supplementary Materials. 


II. STATISTICS OF EIGENFUNCTION 
AMPLITUDES 

As the off-diagonal matrix elements in Eq. Q are 
small, one can employ the perturbation theory for com¬ 
puting the distribution function of the amplitudes x = 
A|^/>(ro)p. The first order perturbation theory gives: 

=\Hnm?/{Hnn- HmmY-, (n ^ to) ( 2 ) 

where the maximum of i/nix) is supposed to be at r = r^. 

The perturbative series converge absolutely if the typ¬ 
ical off-diagonal matrix element |(i7nm)typ| \N~oN 
times the coordination number N is much smaller than 
the typical difference of the diagonal matrix elements 
|(77nn-77mm)typ| ~ W > 5{N). Thus it converges abso¬ 
lutely for 7 > 2 irrespectively of the statistics of diagonal 
matrix elements. For 7 < 2 the convergence of the series 
occurs only because of the random and independently 
fluctuating signs of i7nm and {Hnn — Hmm)- Although 
it is hard to prove such a convergence rigorously, a plau¬ 
sible argument in its favor is that the effective coordi¬ 
nation number of oscillatory contributions is ^fN rather 
than N. The corresponding criterion of convergence is 
A N~oN y/J/ ^ \Y which is satisfied at 7 > 1. 

Consider the regular part of the characteristic func¬ 
tion Q{^) = (e®^). For the Gaussian distribution 
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of matrix elements Eq. Q we obtain Q{^) = Q{^Na), 
where: 

Q(C) = e-</2 Erfc (v'-*C/2) « 1 - \/- 2 </ 7 r. (3) 

The function P{x) = QiC) at a; > 

A^cr ~ 0{N^~'^) is dominated at 7 > 1 by small ('‘Cl. 
That is why it is only the expansion of Eq. (|^ at small C 
what matters for P{x) at 7 > 1. Thus we obtain for the 
regular part of the eigenfunction distribution P{x): 

V{x) = ( v ^ tt )”^ ( 4 ) 

There are two normalization conditions for P(x): the 
normalization of probability Eq. ([^ and the normaliza¬ 
tion of the wave function Eq. (§: 

P{x) dx = 1, (5) 

xP{x) dx = 1. ( 6 ) 

Eq. ^ imposes a cut-off x^^^ ~ to Eq. Q 

at small x, while Eq. Q determines the upper cut-off 
a^max ~ N^~^. A caution, however, should be taken: by 
normalization |' 0 (ri)p = 1 the amplitude |V'(ri)P < 1 
on any lattice site cannot exceed 1 , and therefore x < N . 
One can see that the above estimation for ccmax is valid 
only for 7 < 2 when N'^~^ ^ N. For 7 > 2 a correct 
a^max = N. In order to compensate for the deficiency of 
normalization in Eq. (§ one has to assume a singular 
part of P{x) = V{x) + A5{x — N). One can see that for 
7 > 2 Eq. ^ is dominated by the singular term, and A = 
N~^. This corresponds to the strongly localized wave 
functions. The mechanism of emergence of the singular 
term at the AL transition at 7 = 2 is somewhat similar 
to the Bose-condensation, where the singular term also 
appears because of the deficiency of normalization of the 
Bose-Einstein distribution. 

One can express the distribution function Eq. Q 
through the spectrum of fractal dimension^^^. 



/(a) = lim f{a,N) = lim \n[xNV{x)]/InN, (7) 

N-f-oo N^oo 

where a = 1 — lna;/ln7V or \ijj{ro)\‘^ = N~°‘. Using 
Eqs. Q-0 one obtains: 

f {cf) = 0/2 -t- 1 7/2? (o^min < O < Oniax)- (8) 


The upper cutoff Omax = 7 corresponds to the lower 
cutoff aimin' The lower cutoff amin depends on 7. In 
the localized region 7 > 2, Fig. l]^a), Omin = 0. At 
the AL transition point 7 = 2 the function /(a) has 
the same triangular shape as at lU = Wat on RRG, 
Fig. [3b). In the region of the extended non-ergodic 
states 1 < 7 < 2, Fig. [^c), ccmin = 2 — 7 > 0. It 
is remark able that in the entire region 1 < 7 < 2 the 
symmetrjf^nili] f(^i _|_ 2;) = ^ ^ holds. Finally, at 



FIG. 1: (Color online) The spectrnm of fractal dimensions: 
(a) the singular spectrnm in the localized phase at 7 > 2 . It 
corresponds to the same exponent Tq as for /(a) shown by the 
dashed line, (b) the triangnlar spectrnm at the localization 
transition point 7 = 2 . (c) the spectrum with the gap Omin = 
2 — 7 for the intermediate phase 1 < 7 < 2; (d) the ergodic 
transition at 7 = 1 corresponds to the collapse of Qmax — 

Omin = 2(1 - 7). 


7=1 the two limits Omin and Omax collapse in one point 
a = 1 which marks the transition point to the ergodic 
state. Fig. m- 

Note that f{a) for 7 > 2 (see Fig. [^a)) has a singular 
peak at a = 0 which corresponds to the singular term 
N~^S{x — N) in P{x). This singular f{a) is not a limit of 
any convex function. However, one may easily see that all 
the moments ~ have the same exponents 

Tq as for the “convex” f{a) shown by the dashed line in 
Fig. 3 (a): /(a) = a/7 for 0 < a < 7. Such a triangular 
/(a) with the slope smaller than 1/2 also holds in the 
localized phase on RRGP. 

The numerical calculation of /(a) (see Sec. [^ in the 
Supplementary Materials) which involve the rectification 
and extrapolation procedures described in Ref.®, fully 
confirms the above results. In Fig. [^ we present the 
results of this calculation for N = 10° — 10^"* and the 
extrapolated /(a) (shown by a solid red line) for 7 = 3 
which perfectly coincides with the prediction of our per¬ 
turbative analysis above. The similar coincidence was 
obtained for 7 = 1.5, while for 7 = 0.75 the distribu¬ 
tion function P{x) is practically indistinguishable from 
the Porter-Thomas distribution of the GOE. 


III. THE SUPPORT SET DIMENSION Di 

By calculating the Legendre transform of /(a) 
(shown in Fig. [3 one finds that in the intermediate phase 
1 < 7 < 2 all fractal dimensions Dq = T„/{q — 1) for 
q > 1/2 are the same and equal to (see Sec.jBjin Supple¬ 
mentary Materials for details): 

Dq=D = 2-^, {q>l/2). (9) 

Thus the support set of a typical wave function is a frac¬ 
tal containing sites. As — >■ 00 in the 
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FIG. 2: (Color online) Spectrum of fractal dimension for 7 = 
3 obtained numerically as in Ref.®. The linear part of the 
extrapolated /(a) (solid red line) is exactly as expected /(a) 
(black dashed line). The curves for f{a,N) for increasing N 
are shown by black, blue, green and orange lines from bottom 
to top. The top of the singular peak at a = 0 shown by 
the points of the corresponding color, extrapolates to zero as 
expected (see inset); (inset) the l/\nN extrapolation of the 
singular peak value fpeak = /(O, N). 



FIG. 3: (Color online) Spectrum of fractal dimension in the 
intermediate phase for 7 = 1.5 obtained numerically as in 
Ref.®. All notations are the same as in Fig. [ 2 ] Line colors 
correspond to the same values of N as in FigTbj Expected 
/(a) is shown by a black dashed line. 

limit —)■ 00 it is an extended state. However the sup¬ 
port set contains a fraction of all sites F = tending 

to zero in this limit. Thus it is a non-ergodic state. 

In the localized phase 7 > 2 (including the critical 
point 7 = 2 ) we obtained: 


1.0 



FIG. 4: (Color online) Finite-N spectrum of fractal dimen¬ 
sions f{a,N) for 7 = 0.75 and N = 2^° — 2^® obtained 
by the rectification procedure of Ref.®. For comparison we 
also present f(a,N) (shown by a red line) for the Porter- 
Thomas distribution of wave function amplitudes in the GOE 
obtained by the same procedure at = 2^®. It almost coin¬ 
cides with the (violette) curve for f(a,N) computed at the 
same N = 2^® for our model with 7 = 0.75. In the inset: 
lnP(a;) vs In a; for the same 7 = 0.75 and system sizes as in 
the main plot. The corresponding curve for GOE is shown in 
red. All the curves are almost indistinguishable. 
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EIG. 5: (Color online) The support set dimension Pi (7) and 
the global curvature C'i/iv( 7 ) extracted from the fit (xlna;) = 
(1 — Pi) In A -I- Co -f A“^ vs. 7. The dashed line is the 

prediction for Pi, Eq. (inset) The average ( 2 : In a:) vs In A 

for 7 from 0.75 (bottom) to 2.25 (top) with steps 0.25. It is 
related with the Shannon entropy — IV’(’')I^ In l'*/'(’')P = 

In A — (a; In a;). The global curvature Ci/jv changes sign at 
the transition points 7 = 1 and 7 = 2 . 


-(9) = 


79-1, q< 1/7 
0 , q> 1/7 


( 10 ) 


One can see that the fractal dimensions Dq = 0 only 
at q > 1 / 7 , while they are non-zero and negative for 
0 < <7 < 1 / 7 . This is not exactly the behavior of the 
typical Anderson insulator where all fractal dimensions 
with q > 0 are equal to zero. The behavior similar 
to Eq. (10 1 were found in certain two-dimensional ran¬ 
dom Dirac model^^lHIZl Such a quasi-localized phase is 
referred to as the frozen phase and the corresponding 
transition is known as the freezing transition. In such 
a phase a typical wave function amplitude has several 


sharp peaks separated by valleys where |' 0 p is not expo¬ 
nentially but only power-law small in N (I'f/’ltyp > 
in our case). The same behavior is also found for the 
RRCpwith W > Wat- 

In order to check the existence of the intermedi¬ 
ate phase numerically we computed the average (x In x) 
which is directly related with the Shannon entropy and 
the dimension Di of the support set of fractal wave 
function^^. The results are shown in the inset of Fig. 
where N span from 256 up to 32768. The correspond¬ 
ing values of Di extracted from the linear in In A fit 
are shown in Fig. which are consistent with the tran- 
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FIG. 6 : (Color online) The overlap correlation function 
K{uj) = NJ2r{\''pEir)f (a) For 1 < 7 < 2 func¬ 

tions K{uj)N^ ' for different N collapse into the same curve 
in coordinates Lo/Exh oc lu/N^ ' (see also an inset to Fig.[^. 
(b) For 7 > 2 the collapse of K{u})/N occurs in coordinates 

sitions at 7 = 2 and 7 = 1 . The deviation from the 
expected Di =2 — 7 shown by a dashed line in Fig. 
is a finite-size effect. Indeed, the correlation volume Nc 
close to the localization transition at 7 = 2 is exponen¬ 
tially large Nc oc This follows from Eq. ( [l^ 

of Sec. [V] where the Poisson limit is reached only for 
]V 7-2 ^ 0(1) or In In ~ 1/(7 ~ 2) (see Sec. of 
Supplementary Materials for details). The similar expo¬ 
nential dependence Nc ~ of the correlation 

volume was obtained on the Bethe lattic^^. For system 
sizes N ^ Nc one should see the properties of the crit¬ 
ical point 7 = 2 where D\ = 0. Thus in the vicinity of 
the transition point the support set dimension extracted 
from the finite-size simulations should show a tendency 
towards smaller values as in Fig. However, for 7 < 1.5 
at our systems sizes the support set dimension Di ap¬ 
proaches the values expected in A^ —> 00 limit (dashed 
line in Fig. [^. This fact implies that for 7 < 1.5 we 
reached N ^ Nc and thus it may serve as numerical 
evidence of convergence and existence of non-ergodic ex¬ 
tended phase in the thermodynamic limit. 

We also introduced the 1/Af corrections to the fit with 
its magnitude Ci/j^ being a measure of the global cur¬ 
vature of the (a;lna;) vs. InA^ dependence (see Sec. [C] of 
Supplementary Materials for details). Remarkably, 
changes sign at both the transition points 7 = 1 and 
7 = 2 (though the positive Ci/j^ is very small for 7 > 2 ). 
We also checked that it changes sign at the localization 
transition point of the 3D Anderson model (not shown). 
We believe that the changing of sign of Ci/jy is a conve¬ 
nient way to identify the points of both localization and 
ergodic transitions. 

IV. OVERLAP OF DIFFERENT WAVE 
FUNCTIONS 

Next we compute numerically the overlap of differ¬ 
ent wave functions K{uj) = N J2r{\'^E{r)\‘^\^E+uj{r)\‘^) ■ 


Note that for the ergodic wave functions of GOE K{u!) = 
1 is independent of the energy difference w, as in this 
case the overlap is always 100% (see Sec. of Supple¬ 
mentary Materials for details). Our results presented 
in Eig. 1 ^ show that for 7 > 1 the overlap K^oj) has a 
plateau at w < Eth which is followed by a fast decrease 
K{oj) oc l/w^ for oj > Etk- The Thouless energy Eti^ 
that separates the GOE-like behavior (plateau) from the 
system specific behavior {K{oj) oc W^), depends on N 
as a power-law N~^. However, the scaling exponent z 
is different in all the three phases (see Eig. |^. In the 
localized phase Fig.j^b) we obtained z = 7/2 which cor¬ 
responds to rare resonances when w < ~ N~'^/‘^. 

In the extended non-ergodic phase Fig. ^a) we found 
Erph ~ N^~"^ ~ 6 N^, where S = l/{Np{^ is the mean 
level spacing. This corresponds to all N^ sites in the 
support set being in resonance with each other. The cor¬ 
responding N^ states produced by linear combinations 
of basis states localized on resonant sites form a mini¬ 
band of levels of the width Exh ~ ^ N^. Clearly, the 
states inside such a mini-band should have the GOE-like 
correlations. On the contrary, the states separated by the 
energy distance w > Exh should belong to different sup¬ 
port sets which poorly overlap with each other. At the 
ergodic transition at 7 = 1 and in the entire extended er¬ 
godic state at 7 < 1 we obtain Exh ~ 0{N^) (see Sec.[D| 
of Supplementary Materials for details), and the plateau 
extends to entire spectral band-width. The emergence 
of such a plateau that survives the limit V —>■ 00 is a 
signature of the ergodic stat^^. 

Surprisingly, the overlap function K{uj) ~ N^~'^/ui^ —>■ 
0 as V —>■ 00 at any fixed u and 7 > 1. This phenomenon 
of “repulsion of wave functions”^ is a peculiar feature of 
our model. The non-ergodic fractal states at the localiza¬ 
tion transition points of the two and three- dimensional 
Anderson models of the Dyson symmetry classes, as well 
as those of the power-law banded random matrice^^^MIini 
show a different behavior. In these models, the Thouless 
energy for fractal states is proportional to the mean level 
spacing Exh ~ ^ and the behavior for w > Exh(N) is 
described by the conventional Chalker’s scalin^^iE^Illl 

K{uj) ^ (II) 

Only at a very large energy separations uj of the order 
of the total spectral band-width, the “repulsion of wave 
functions” was observecP^. 

A remarkable feature of the present model is that the 
Thouless energy in the region of extended non-ergodic 
states is much larger than the mean level spacing: 

ExhS N°^ (X N-\ ( 12 ) 

One can interpret this relationship as a non-trivial dy¬ 
namical scaling exponent 

z = 1-D2<1. (13) 

Eor non-interacting systems in two or three-dimensions 
in the point of Anderson transition z = 1 for all Dyson 
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symmetry classes. A non-trivial z is known only in 
two-dimensional systems described by the Dirac equa¬ 
tion with random vector-potential which belong to chiral 
symmetry classe^^lHlZl the freezing transition is ob¬ 
served. 

In terms of the dynamical exponent z the leading 
power-law term in the Chalker’s scaling for w 3> E^h 
can be rewritten a^ 


K{uj)(XU! ^JL = {I - D 2 )/z. (14) 

In our model we have: 

Eth 


K(u 


9 ’ 


(w » Erh)- 


(15) 


One can consider Eq. (15) as a particular case of ex¬ 
pansion in Exh/^ 1 with the leading term exponent 

/i = 1: 


K{uj) 


1 

ui 


Co -I- Cl 



LO 


(16) 


in which the coefficient cq is zero. We will see below that 
Eqs. 16) hold for the RRG too. However in this case 
the coefficient Cg 1 is non-zero. Thus one can speak of 
the special universality class of models with a non-trivial 
z ^ 1 and n = I which the present model belongs to 
together with the RRG model. 


V. SPECTRAL FORM-FACTOR 

Finally, we present the results of a rigorous calculation 
of the spectral form-factor C{t,f ) = e**^"*"*"** 

with a set of eigenvalues {En} of H. To this end 
we generalize to the case 7 > 1 the results of Ref.l^ 
where was derived for URP model Eq. ([^ with 

7 = 2 using the Itzykson - Zuber formula of integra¬ 
tion over unitary group. The final result (see details 
of the derivation in Sec. 1^ of Supplementary Materials) 
for = 2TT5{t + — t')) — 1 ] is given by 

Eq. ([^ 


S{u) = 1 -k ^ 


'2/i(ku3/2) 


KU 


. 3/2 


47r 


X dx 

\Jx + 1 




(17) 


with the modified Bessel function Ii{x), Exh = 
S K = -s/STriW^A^ and A = Ap(0). The unfolded 



FIG. 7: (Color online) Unfolded spectral form-factor R{t) for 
the RP model for two different cases: (i) 7 = 2 and (ii) 7 = 1.5 
at A = 10^, 10®, 10^. The falling part corresponds to attrac¬ 
tion of levels while the rising part corresponds to repulsion of 
levels. The GUE form-factor is shown by the red dashed line. 
The Poisson distribution corresponds to R{t) = 1. (inset) 
The spectral form-factor in the variable u = |(t — t') Eth 
with Exh = 5 A® *' in the A —>■ 00 limit for different val¬ 
ues of 7. There are five distinctly different phases: insulator 
7 > 2, AT critical 7 = 2 , non-ergodic extended 1 < 7 < 2, 
FT critical 7 = 1 and ergodic 7 < 1. 

spectral form-factor R{t) with i = \{t — t')S is given by 
i?(t) = 5'(tA^“'''). It follows from Eq. ( [TtI that for 7 > 2, 
R{i) —>■ 1 in the A —>■ 00 limit, which corresponds to com¬ 
pletely uncorrelated energy levels and the exact Poisson 


statistics. Another important feature of Eq. ( [l7| is that 
R{0) = 1 for all 7 > 1. 

In Fig. [ 7 ] we plot the unfolded spectral form-factor 
R{t) for the two phases: (i) the critical phase of the AL 
transition at 7 = 2 and (ii) the intermediate phase at 
1 < 7 < 2. One can see that while at 7 = 2 the func¬ 
tion i?(t) has a non-trivial A —>■ 00 limit, for 1 < 7 < 2 
the limit coincides with that of the GUE, except for the 
point t = 0 where there is a jump in R{t). This jump is a 
hallmark of the intermediate phase. To demonstrate this 
more clearly we blow up the region of small t by re-scaling 
the variable t =k u = ^{t — t')5 — t') Exh- 

Note that Exh = appears again as the char¬ 

acteristic scale where the level repulsion is taken over 
by the level attraction in Fig. In this new variable 
R{uN^~‘^) = S{u) has a non-GUE A —>• 00 limit: 

S{u) = exp(—27rA^ u), 1 < 7 < 2, (18) 


which is shown in the inset to Fig. The true GUE 
form factor is just identically zero in this limit. The 
existence of the new scale Exh and a non-GUE A —>■ 00 
limit Eq. (18) in the variable u= — Exh have been 
overlooked in Ref.^Sl. 

Eq. (18) holds for u > Exh ~ A®“^, and S{u) is satu¬ 
rated at S'(O) « ^ ^ for u < A^”'’', which corre¬ 

sponds to |t — smaller than the inverse total spectral 
band-width (see Sec. of Supplementary Materials for 
details). For 7 = 1 we have A^“'>' = 1 . Thus the value 
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FIG. 8: (Color online) The support set dimension D\ and the 
global curvature Ci/jv as the function of disorder strength 
W extracted for the RRG model as in Fig. No ergodic 
transition is detected by changing of sign of the global curva¬ 
ture. This is consistent with the statemeniP that the entire 
extended phase is non-ergodic. 


S(0) ^ at 7 = 1 is smaller than 1 . So, in ad¬ 

dition to the specific critical behavior of S{u) at 7 = 2 
(shown by the red curve in Fig. one obtains yet an¬ 
other critical behavior of S{u) at the ergodic transition 
7 = 1 (shown by the dashed yellow line in Fig. which 
is stable in the N ^ 00 limit and is characterized by 
5(0) < 1. For 7 < 1, is increasing with N, making 

5(0) = max5(u) —>■ 0 as > 00 . This is how the GUE 
limit S{u) = 0 is reached. 

Note that the fact that Eth ^ 5 affects the level 
number variance var(n) (n and var(n) are the average 
number of levels and the level number variance in a cer¬ 
tain spectral window, respectively) in which a new scale 
Erh/S ^ appears for n at 1 < 7 < 2 : 


var(n) 


~lnn, 1 <C n <C 
n, N‘^~^ n N. 


(19) 


The level compressibility is ill-defined in our model 

because of the jump in the spectral form factor i?(t)|jv_,.oo 
at t = 0. Formally it can take any value from % = 0 to 
X = 1 depending on the parameter n/N‘^~^. This is 
in contrast to other random matrix models with multi¬ 
fractal eigenstates, e.g. PLBRM and Moshe-Neuberger- 
Shapiro models (see Ref.l^ and references therein) where 
the level compressibility is well defined and takes a defi¬ 
nite value 0 < X < 1 in the N ^ 00 limit. 


VI. COMPARISON WITH RRG MODEL 

In the Introduction we mentioned a heuristic relation 
between the Anderson model on a hierarchical RRG and 
the RP model which has no apparent hierarchical struc¬ 
ture. It is instructive now to compare the main results 
of this paper with the corresponding results for RRG. 

First of all we recall (see Fig. and the corresponding 
explanations in the text) that all the moments (|V’P‘^) in 


FIG. 9: (Golor online) Gollapse of data for K{uj) for the An¬ 
derson model on RRG with the branching number K = 2 
and N = 2000, 4000, 8000, 16000, 32000. The collapse oc¬ 
curs in the coordinates K{uj) and lo /which 
corresponds to Eqs. (121, and Fig. j^a) in RP-ensemble. 
The values of fractal dimension D 2 extracted from the best 
collapse are in excellent agreement with the corresponding 
values found from the moments (see inset). Different 

disorder strengths W correspond to different collapse curves 
and different dynamical exponents « in Eth = N~^. This ex¬ 
cludes a possibility for an apparent non-ergodicity of extended 
phase to be a finite-size effect reflecting properties of only one 
single critical disorder strength. The behavi or of K{uj) ~ l/oj 
for oj > Eth = 5 corresponds to Eq. (161 with a non¬ 
zero Co ~ 1. (inset) The fractal dimensions D 2 = 1 — 2 as a 
function of W found from the best collapse of data for K{lj) 
(red circle) and from the scaling of moments dV'l"^) with N 
(black circle). 


the localized and the AT critical phases of our model have 
exactly the same q-dependence as in the corresponding 
phases of the RRG. The N- and 7 - dependence of the 
moments is also very similar (cf. Fig. with Fig. to 
the corresponding N and VF-dependences for the random 
wave functions obtained by the exact diagonalization of 
the Anderson model on RRG with the branching num¬ 
ber K = 2 and N = 2000 — 16000. However, there is an 
important difference: we found only one point of chang¬ 
ing the sign of Ci/jy on RRG which corresponds to the 
known point of the Anderson localization transition at 
W « 17.5. 

In Fig. [^we demonstrate that in the case of RRG the 
scaling of the Thouless energy with the system size fol¬ 
lows the same Eq. (12) as for the present RP model and is 
different from the standard Chalker’s scaling. The falling 
part of E(uj) at uj ^ Eth can be described by the unified 
expansion Eq. (16) both for RRG and our model, albeit 
the coefficient cq is zero for the present model and is non¬ 
zero for RRG. It is important that for RRG the data 
for properly re-scaled F(y) = N~^K{yN~^) at differ¬ 
ent N collapse on one same scaling curve F{y) = Fw{y) 
that, however, depends on W, as well as the exponent 
2 = 1 — D 2 (see inset of Fig. [^. This is very different 
from the usual scaling E(u!) = N°- F{ujN^, Nc/N) in the 
vicinity of a single critical point W = Wc where the ex- 






















ponents a and z are determined by the property of this 
critical point and not by the distance \W — Wc\ from 
this point which determines only the correlation volume 
Nc{\W — Wc\)- In our opinion, this implies that there 
is a line of critical points at < Wc which determine 
the behavior of the system at least in a parametrically 
large interval of sizes Nc 2 N ^ with the second 
characteristic size scale Nc 2 Nc. 

Our conclusion is that the localized and AT critical 
states are very similar for our model and RRG. The ex¬ 
tended states show non-ergodicity in a broad interval of 
7 and W and are characterized by the Touless energy 
which in both models is much larger than the mean level 
spacing. However, the existence or non-existence of the 
ergodic transition is more subtle and depends on tiny fea¬ 
tures of the model. It exists in our model and most prob¬ 
ably does not exist on RRG with the branching number 
K = 2. Nonetheless our study largely confirms expecta¬ 
tion on the similarity between the RRG and RP models. 
This allows us to speak on the special class of models with 


the explicit (RRG) or hidden (RP) hierarchical structure. 
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Appendix A: Numerical extrapolation of the spectrum of fractal dimensions. 


We start with the computations of the distribution function P{x) of the normalized amplitude x = N\il}\’^ in the 
Rosenzweig-Porter (RP) ensemble and extracted the spectrum of fractal dimensions f{a,N) = ln(A^Pe„^(lnA^(l — 
Q;)))/lniV, using the approach of Ref.® In order to eliminate the effect of zeros of wave functions ip which dominate 
the distribution function P{x) at small x < Xmin ~ N^~'^ and extract the distribution function of a smooth envelope 
of |'0P we represent ip = ipenv x Vi where r] is the Gaussian Orthogonal Ensemble (GOE) random oscillations with 
the unit average square which are supposed to be statistically independent of ipenv ■ Then the distribution of In a; is a 
convolution of the distribution Penviv) oi y = \nxenv = and the known GOE distribution of In 77 ^. Making 

a numerical de-convolution one obtains the distribution Penviv) in which the effect of zeros of ri{ro) is eliminated. 
Such a “rectified” distribution function decreases much faster at small Xenv < Xmin than the distribution P(lna:). This 
“rectification procedure” results in a sharp cut-off of /(a, N) at large a. Then the so obtained f{a, N) is extrapolated 
to TV = 00 using the ansatd^ f{a,N) = f{a) + c{a) / In N. 

The results for the localized case 7 = 8 and N = 2®...2^^ are shown in Fig. One can see that f{a) rectified and 
extrapolated as explained above has a linear in a part which exactly coincides with the prediction of the paper (see 
Fig. 1(a) of the paper). Moreover, the singular peak predicted for 7 > 2 in the paper is observed and its top value 
tends to zero as N ^ 00 (see inset in Fig. [^. In Fig. we present the extrapolated /(a) for 7 = 1.5. Again, the 
linear part of f{a) coincides with the expectation shown in Fig. 1(c) of the paper. 

Finally, we present the results for f{a,N) for 7 = 0.75 < 1 (see Fig. [^. One can see that the distribution function 
P{x) is weakly iV-dependent and is almost indistinguishable on the logarithmic scale from the GOE distribution 
function. The function f(a,N) is more sensitive to N but also in this case there is a minor difference from the 
corresponding GOE result. 


Appendix B: Moments of \ip\‘^. 


Here we discuss the moments 


I,= 



N{\iP\^‘^) cx 


(Bl) 


We start by considering the localized phase 7 > 2 where the moments with q > 1 are dominated by the singular 
term N~^S{x — N) in the distribution function. One can easily see that all such moments are A-independent, which 
corresponds to T{q) = 0. Using this equality and given that T{q) and /(a) are related by the Legendre transforirP^l 

T{q) = qaq - /(a,), /'(a,) = 9, (B2) 


and that the singular peak is located at a = 0 , one immediately obtains tha t the top of the peak corresponds to 
/(O) = 0. Considering also 9 < 1 and using the Legendre transform Eq. (B2) we obtain from /(a) of Fig. 1 of the 
paper: 


T-{q) 


( 79 - 1 , 9 < 1/7 

10, 9 > 1/7 ’ 


(7 > 2). 


(B3) 


By the nature of the Legendre transform the tangential to /(a) with any slope 9 is the same for the singular and 
“non-convex” /(a) shown by the red solid line in Fig. 1(a) of the paper and for the “convex” f{a) shown by the blue 
dashed line. Thus for both types of f{a) shown in Fig. 1(a) of the paper T{q) is given by Eq. (B3). 

For 1 < 7 < 2 Eq. (B2) gives: 


^( 9 ) {(2-^7)(g-l), g> 1/2’ (^ < O'< 2)- 


(B4) 


Eq. (B4| implies that in our model the multifractal dimensions Dq for 1 < A < 2 do not depend on 9 for 9 > 1/2 and 
are equal to: 


Dq = T{q)/{q-l) = (9 > 1/2). 


(B5) 
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Appendix C: Support set dimension and the curvature. 

We computed numerically the moment {x Inx), where a; = iV which in a pure multifractal state should behave 
as: 


{x Inx) = (1 — Hi) InTV + const, {N —>■ c»), 


(Cl) 


where Di is a dimension of the wavefunction support set. In an ergodic phase Di = 1, while in the localized phase 
Di — 0. In the intermediate multifractal state 0 < Hi < 1. The support set dimensiorP^ may be expressed through 
the solution a = oi of the equation: 

f\a) = 1, (C2) 


(with /(a) being a spectrum of fractal dimensions) i.e. Hi = oi is a point a = ai where the line with the slope 1 is 
tangential to f{a). One can immediately calculate Hi as a function of 7 using Fig. 1 of the paper: 


1, (7 < 1) 

Hi = <(2- 7 , ( 1 < 7 < 2 ) 

0 , (7 > 2) 


(C3) 


The dependence of the moments Ii = {x Inx) vs. InIV is shown in the inset of Fig. in the main text for 
N = 2®, 2®..., 2^®. By fitting to Eq. (Cl) one can find Hi and compare it with the expected result Eq. (C3). The 
comparison is shown in Fig. in the main te xt. One can see that the apparent Hi extracted from limited sizes 
A = 2® — 2^® deviates from the prediction Eq. (C3) as 7 approaches the Anderson transition point 7 = 2 . We believe 
that this is a finite-size effect related to the correlation volume W ~ that diverges exponentially at the transition. 
Clearly, for N < A)., one should see the properties of the critical point where Hi = 0. That is why the apparent Hi is 


smaller as the prediction in the vicinity of 7 = 2. At the same time, Eq. (C3) describes very well the data points for 


Hi close to 7 = 1. This probably means that the ergodic transition is not associated with an exponentially divergent 
correlation length 


One can quantify the deviations from the linear behavior Eq. (Cl) by introducing the 1/N correction which has a 
finite curvature in the In A variable. In reality, the finite-size scaling exponents are not known for this model, and 
the true corrections could be very different from 1/N. However, coefficient Ci/n in the simple fit: 


(x Inx) = (1 — Hi) In A-I-Co-I-Ci/iv^ (C4) 

gives an idea about the “global curvature” of the dependence (x Inx) vs. IniV. The dependence of this coefficient 
on 7 is shown in Eig. [^in the main text. It has a characteristic peak shape. For 7 in the vicinity of 1 the finite-size 
effects are small at our system sizes, and the global curvature is small too. Close to the point 7 = 2 our system sizes 
are far too small to deviate from the critical behavior at the transition point. In this case the global curvature is 
small too. The absolute value of the global curvature reaches its maximum where In A ~ In A^c- 

An important observation is that the “global curvature” changes sign at the ergodic transition. This observation 
may be used to locate the transition point in finite-size calculations. 

In addition to the global curvature one may also introduce a “local curvatures” as follows: 


(Ji[[l]]+/i[[3]]-2/i[[2]])/ln^2 (/J[_i]]+j,[[_3]]-2Ji[[-2]])/ln^2 

(1 + {h[[l]] - /i[[3]])V41n^ ^ (1 + m-l]] - /i[[-3]])V41n2 2f^ ’ 


where is a moment Ii for the i-th system size counted from the firs t one N = 2®, while /i[[—i]] is the moment 

/i for the Ath system size counted from the last one N = 2^®. Eq. (C5| is nothing but the discrete variant of the 
curvature: 


K = 


h"{x) 

(1 + [H(x)]2)3/2 


(C 6 ) 


of a curve given by a function h{x). 

In the table below we present «:< and for different 7 . One can see that the local, as well as the global curvature 
is negative for 1 < 7 < 2. It is positive and small for 7 = 0.75 < 1. Thus it is likely, that the local curvature, 
too, changes sign close to the ergodic transition at 7 = 1. The absolute value of k decrease with increasing N in the 
vicinity of 7 = 1 (for 7 = 1.25 and 7 = 1.40), which signals about convergence. In contrast |k| increase with increasing 
N in the vicinity of 7 = 2 (for 7 = 1.6 and 7 = 1.75) as the system size N approaches the correlation volume from 
below. For 7 = 2.0 and 7 = 2.25 the local curvature is very small and is inside the error bar. 
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TABLE I: Local curvatures k< and k>. 



ft< X 10 

K> X 10 

7 = 0.75 

-f0.05 

-h0.05 

7 = 1.00 

-0.16 

0.03 

7 = 1.25 

-0.42 

-0.10 

7 = 1.40 

-0.40 

-0.27 

7 = 1.60 

-0.13 

-0.17 

7 = 1.75 

-0.03 

-0.06 

7 = 2.00 

0.00 

0.00 

7 = 2.25 

0.00 

0.00 


Appendix D: Overlap correlation function and the Thouless energy. 


Conventional multifractal correlations imply not only a power-law scaling of the moments of |'(/'ra(fo)P with the 
system size N but also a specihc power-law two- and multi- point correlations. In particular, the overlap correlation 
function: 


K{uj) = N'^\iljE{r)\^\i}E+u^{r )\^, /x = (1 - £> 2 ), 


(Dl) 


for 6 < oj < Eq (wit h <5 bein g the mean level spacing an d E q ~ 0(1) being the onset of the anti-correlations) obeys 
the Chalker’s scalin^^SIlIIlll energy domain. Eq. (Dl) holds exactly at the critical point of the 3D Ander son 

model and in certain random matrix ensembles, e.g. for the power law banded random matrices (PLBRM)p2inil. it 
implies an enhancement of correlations compared to the case of independently fluctuating wave function which would 
result in K{uj) = 1. However, Eq. (Dl) is also approximately valid in the metallic phase of the 3D Anderson model 
close to the transition pointPS in the limited range of dj < uj < Eq, where = {pNc)~^ is the mean level spacing in 
the correlation volume Nc ~ For to < the correlation function saturates developing a plateau, and for uj > Eq 
it decreases as fast as 0 ;“^^. In the 3D Anderson model the plateau survives the thermodynamic limit N ^ 00 and 
extends to larger w as one goes deeply into the metallic phase. It is thus a signature of the ergodic extended state. 
In the region lo > Eq, the overlap correlation function is small, which signals on the “repulsion of wave functions” at 
large energy separations^. A similar phenomenon of eigenfunction repulsion for uj > Eq was observed in the PLBRM 



FIG. 10: (Color online) Overlap correlation function Eq. (Dl I: (a) for the Rosenzweig-Porter model at different 7 at A = 2^“^; 
(b) for the PLBRM model (red symbols),GOE (green symbols) and banded random matrix model (blue symbols). There are 
three regions with different behavior for the PLBRM: the plateau, the Chalker’s scaling with the non-trivial exponent , 

the fast decay ~ Only two of these three regimes are present on the plot (a): the Chalker’s scaling is absent. The GOE 
behavior K{u)) = const is identical to the one for RP model at 7 = 0.75. Note that the behavior in the localized region 7 > 2 
of the RP model is qualitatively different from that of the quasi-one dimensional localization in the banded random matrices: 
in the former case the repulsion of wave functions is present, while in the latter case positions of centers of localization are 
randomly distributed in space with almost no correlations. 
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with a small bandwidth b (see^or Fig. 10’b)). 

We calculated K{u}) in our model numerically. The result is presented in Fig. 3 of the paper and in Fig. [Io|[a) in 
the SM. Surprisingly, no Chalker’s scaling and enhancement of correlations similar to Eq. (Dll was observed. For 
1 < 7 < 2 the correlations are fast decreasing at a; > Exh ^ 


K{uj) 


uj>ETh 

N^-\ u<ETh' 


(1 < 7 < 2) 


(D2) 


like in the high-energy region a; > Eq of the 3D Anderson and PLBRM models. The GOE-like plateau is present 
only in a narrow interval of small w < Eth which shrinks to zero in the thermodynamic limit. Thus we 

can identify Eth N'^~^ as the Thouless energy for this model, which by definition is the border line between the 
GOE-like behavior (plateau) for uj < Exh and the system-specific behavior for uj > Exh (fast decay oc W^). 

Note that the onset of the anti-correlations scales with N exactly as the parameter Exh = which enters the 

dimensionless variable u = — Exh in the unfolded spectral form-factor S{u) in the level statistics (shown in the 

inset of Fig. [^in the main text). It can be nicely expressed in terms of the number of populated sites in the 
wave function support set: 

Exh = SN^\ (1<7<2). (D3) 

Qualitatively similar (but quantitatively different) behavior is observed in the localized phase 7 > 2: 


In this case we have: 


K{uj) 


uj> Exh'^ N-^/^ 

N, UJ < Exh 


(7 > 2 ) 


(D4) 


Exh ~ (7 > 2). (D5) 

Surprisingly at 7 > 1 and a fixed w the overlap K{uj) drops below the independent wave function limit K{uj) = 1. 
Thus in N ^ 00 limit of our model for 7 > 1 the repulsion of wave functions happens at all energy scales. 

For 7 < 1 the plateau ^ 1 in K{uj) extends to almost the entire spectral bandwidth (see Fig. [To|(a)) and may even 
exceed it: 


Exh ~ 0(A^°), (7 < 1 ) 


(D 6 ) 


The emergence at 7 = 1 of a plateau that survives the thermodynamic limit and occupies a finite fraction of (or all 
of) the spectral bandwidth is a very clear signature of the ergodic transition at 7 = 1 . 


Eqs. (D2 D4) were checked by a collapse of the data points for a fixed 7 but different N = 2^^ — 2^^. The results 


are shown in Fig. [^in the main text and Fig. EH „ 

As 7 approaches the Anderson transition point 7 = 2 deviations from Eq. (D2| emerge. Eor 1 < 7 < 2 the best 
collapse was found to occur in the coordinates X = Y = K{uj)/N^^^'^>, where: 


Xx{n)=l-l + d, Ay( 7 ) = 7 - 1 - 2 d. 


(D7) 



FIG. 11: (Color online) Collapse of data for K(uj) at 7 = 2.00 and 7 = 2.25 in the coordinates X = uj/N Y = K{u])/N. 
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FIG. 12: (Color online) Collapse of data for K{lo) at 7 = 1.25 and 7 = 1.50 in the coordinates X = uj/N^ , Y = 
K{lo)N^ corrections d = 0.04 and d = 0.10 are obtained from the best collapse. 



FIG. 13: (Color online) (a) The support set dimension Di vs. 
finite-size effects make Di and 1 4 - Ax closer to its value at 7 = 



7 . (b) The scaling exponent Ax (7) vs. 7 . In both cases the 
: 2 than expected (dashed line). 


The correction d{j) is most probably due to the finite-size N that could be comparable with the correlation volume 
near the Anderson transition point 7 = 2 . The dependence 1 -I- Ax( 7 ) vs. 7 is plotted in Fig. [I^b). It has an 


apparent analogy with the dependence Hi( 7 ) of the support set dimension in Fig. 13’a). We believe that in both 


cases this is a finite-size effect when N < gcorist/fz 7) jg gj^aller than the exponentially large correlation radius near 
the AT point 7 = 2 . 


Appendix E: The spectral form-factor. 


Here we consider the spectral form-factor ({H„} is a set of eigenvalues of H) which 

was derived for URP model Eq. (1) in the main text using the Itzikson-Zuber formula of integration over unitary 
group and we start by Eq. (3.2) of Ref.E^ 




Ntt' 


dz 

2TTi 


dz' 

2TTi 


z') 


9{z,z') 


/\N 


1 + 


rr 


{z' — z — t){z' — z t') 


'\N 


-p{zrp{z') 


■ (El) 


In this exact formula integration is extended over the contour Fjj that encompasses the real axis and r = ita. The 
functions g{z^z') and p{z) depend on the diagonal disorder distribution p{a) through a{z) = J p{a)da/{z — a) = 
{{z — a)~^) and are defined as follows: p{z) = 1-1- Ta{z), 


giz,z') 


1 -I- Ta{z) + T'a{z') 



[a{z) - a{z')]. 


(E2) 


In order to do the limit of infinite matrix size A —> 00 we make a re-scaling: 


{z,z') 


- ± s 
2 

-L y 


ie 




(E3) 

(E4) 
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where e —>■ +0 and q,q' = ±1 for the part of the contour below and above the real axis. The exponents At > 0 
and > 0 should be chosen so that (i) finite limits limjv-too[ 5 ( 2 :, and limjv-too[p(- 2 )]'^ exist, and (ii) the entire 
expression is finite in the iV —?> 00 limit. One can show that the choice 

At = 7 - 1, Az = 1, (E5) 


satisfies all these conditions if 7 > 1. Now doing the limit A —>■ 00 at fixed T, s, x, y we observe that the parameter 7 en¬ 


ters the limiting expression Eq. (El) only in the exponents e 

2 


gz(tz-l-t'z') ^ ^iTx+iN-^- 


This observation immediately tells us that for 7 > 2 the spectral form-factor C(t,t') vanishes, which leads to 
completely uncorrelated energy levels and the exact Poisson statistics. In the case 7 = 2 considered in Ref.^ the 
level statistics is different from Poisson, as C{t,f) = 2ttp{0)6{T) [S'(m = s/p{0)) — 1] is not zero. However, the level 
compressibility equals unity 


X = ^(0) = 1. 


(E 6 ) 


like for uncorrelated energy levels. 

One can easily see that Eq. (E 6 ) remains valid also in the entire region 1 < 7 < 2, since at s = 0 the 7 -dependent 
expone nts e~^ ® and are equal to 1 anyway. 

Eq. (El) can be cast as follows + K 2 {t,t')), where: 




1 


Ntt' 




2Tri 


dz 


2'Ki 


{g{z,z'Y - p{zYp{z')^) , 


(E7) 


1 


dz 


= — (t> (p ::^e 


dz' 


9{z,z') 


i\N 


z') 

N ./p^ 27ri Jp^ 27rj iz' — z — t){z' — z + t’') 


(E 8 ) 


Performing the re-scaling Eqs. (E3 - |E5[ ) of the paper and changing the variables y ^ {q — q')y we arrive at: 

-11 , (E9) 




dx 


K2it,t')=J2 / .— / —e 


dx 


ixT—2-KX^sqp{x) f ^ib(y—ie) _ y 

27r J 47 r {y — ie + isqX^/2)'^ 


-S If 


ixT 


<^y,ibv 1 

47r {y + isqX^l2Y ’ 


(ElO) 

where p(x) is the distribution of Hnn coinciding with the density of states for 7 > the limit A —>■ 00 is taken 
everywhere, except in the exponents e*^*^"*"* ^ ^ —>■ ^ summation over g' is taken, and: 


a = 7 rA"^s"^p(a:), b = 2sqN~' 


-2 


(Ell) 


The exponents determine the allowed contour deformations. Deforming the contour so that the exponents are small 
at \y\ —0 we observe that both Ki{t,t') and K 2 {t,t') are identically zero for sq < 0. Thus the summation over 
q can be dropped with s being replaced by |s|. One can also express the integrals over y in terms of the Bessel 
functions, deforming the contour so that it encompasses the origin along the infinitesimal circle. The final result for 
C{t, t') = 2TT5(t -|- t') [5'(u = s/p(0)) — 1] reads: 


S{u) = 1 + 




2 e- 






(E12) 


47r 


: dx 


'Jx + l 


{x + 1 )) e 


-xu^A'^N'’-^ 


coinciding with Eq. (9) in the main text. Here Ii{x) is the modified Bessel function and A = Xp(0). Eq. (E12) is valid 


for 7 > 1. The second and the third terms in Eq. (E12) correspond to Ki{t,t') and K 2 {t,t'), respectively. At 7 > 2 


both of the terms vanish in the A —>■ 00 limit, and the Poisson limit S{u) = 1 is reached; for 7 = 2 both the terms 
are non-zero; for 1 < 7 < 2 the third term cancels the first one and the second term tends to a finite limit as A —>■ 00 : 


lim S (u = l(t — t')S '*') = e 

N^oo V 2 \ / 


— 27rA^ u 


(E13) 
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The unfolded spectrum form factor is given by: 

i? (t = - t')6) =S{u = , 


(E14) 


where 6 is the mean level spacing. 

For 7 > 2 one can see that R{t) tends to the Poisson limit R{t) = 1 a.s N —>■ oo. 

For 1 < 7 < 2 and any t > 0 it evolves towards the Gaussian Unitary Ensemble (GUE) form factor: 


R, 


GUE 


(t) = 


t/(27r), 0 < t < 27r 
1, t>27r 


(E15) 


as N ^ oo. 

However, in contrast to RcuEi^) the function i?(t) has a jump at t = 0: 

R{t = 0) = 1. 


(E16) 


On the other hand, the function S{u) = R{t = uN"^ at 1 < 7 < 2 has a finite non-singular limit Eq. (E13) as 
N —>■ 00 : 


lim RluN^-^) = e 

N^oa 


— 2\ 27rA"^ n 


while for the true GUE form-factor this limit is zero: 


lim Rque{uN^ ) = 0 . 
N—¥00 


The limit Eqs. ( E13|E17 ) and the new emergent energy scale (the Thouless energy) 

Erh = 5N^-^ ~ N^-^. 


(E17) 


(E18) 


(E19) 


that separates repulsion of energy levels at small energy difference (large “time” difference t) from attraction of energy 
levels at large energy difference (small t), is a hallmark of the extended non-ergodic state for 1 < 7 < 2 . 

For 7 < 1 Eq. (E12) does not hold. The formal reason is that in the re-scaling Eq. (E3) the second term is 
no longer the leading one as N —>■ 00 . The physical reason is that at 7 < 1 the Thouless energy Eq. (E19) is no 


longer small compared to the spectral band-width EdoS (the width of the mean density of states p{e)). Note that 
S{u = ^{t—t')E'rh)j which is the Fourier transform of the two-level correlation function, must saturate when {t—t')ETh 
gets smaller than Exh/EuoS^ where EjooS is the maximal energy difference w within the spectral band-width: 


S{u) = e 


_ „ — ETh/Eoa 


{u <C Exh!E dos)- 


(E20) 


If Exh/Ex)oS is formally divergent, as it is the case for 7 < 1, we conclude from Eq. (E20) that Smax —t 0. This is 
how the true GOE limit S{u) = 0 is reached at 7 < 1. 

Note that the fact that Exh ^ ^ affects the level number variance var(n) (n and var(n) are the average number of 
levels and the level number variance in a certain spectral window, respectively) in which a new scale Exh/S ~ N‘^~^ 
appears for n at 1 < 7 < 2 : 


^ Inn, 1 <C n <C 
n, n<S^ N. 


(E21) 


The level compressibility x = lim]v->cxj, = 0 if limw^oo, n/7V^ ^ = q and y = 1 if limw^oo, n/N'^ '>' = 00 but 

limw-nxj, n/N = 0. 










